hep-ph/9406252 

KEK-TH-376 

KEK preprint 93-160 

UT-659 

March 1994 
H 



Gluon and Charm Distributions in the Photon 

K. Hagiwara, M. Tanaka and I. WatanabeB 

Theory Group, KEK, Tsukuba, Ibaraki 305, Japan 

T. Izubuchi 

Dept. of Physics, Univ. of Tokyo, Bunkyo-ku, Tokyo 113, Japan 

(March 1994) 

Abstract 

We study systematicahy the sensitivity of the photon structure function data 
on the gluon contents of the photon, by using the leading order Altarelli-Parisi 
(AP) equations. Charm quark contribution is studied in the quark parton 
model and by using the massive quark AP equations of Gliick, Hoffmann and 
Reya. The present photon structure function data are found to allow wide 
range of gluon distributions in the photon. We give a set of six scale-dependent 
parton distributions in the photon (WHITl to WHIT6), which have system- 
atically different gluon contents. Sensitivity of the structure function at small 
X and that of the total charm quark production cross section to the effective 
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gluon distribution are discussed. 
12.38.-t 
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I. INTRODUCTION 



The deep structure of the photon has been expected to be calculable perturbatively in 
QCD [jl|, except at small x [Q. In practice, however, non-perturbative effects are found to be 
significant at present experiments, where the electromagnetic structure of Weizsacker- 
Williams quasi-real photon ^ has been measured up to the momentum transfer scale 
= 100 GeV^ in e^e~ collision experiments. Several parametrizations of the scale-dependent 
effective parton distributions in the photon have been proposed; some are based on plau- 
sible dynamical assumptions and the others P|,pi| have been obtained by fitting 
phenomenologically to the photon structure function data [|Tl1-[l7||. These parametrizations 
typically have similar quark distributions which are directly constrained by the structure 
function data. On the other hand, wildly different gluon distributions have been proposed 
since the electromagnetic structure of the photon is rather insensitive to its gluon content. 

Recently TRISTAN experiments [p!8| , p!9| have shown evidence for the effective gluon con- 
tent of the photon in two-photon production of high transverse momentum (px) jets. The 
observed jet production cannot be understood without the contribution from gluons in the 
colliding photons, whereas it does not allow a very hard gluon distribution [|1^] that is 
consistent with the structure function data. 

More recently, the TRISTAN experiments have reported evidence for copious production 
of charmed particles in the two-photon collision process, in their various decay modes and 
at various pt range pO|-p^. The charm production rate has been found to be particularly 
sensitive to the gluon distribution in the photon ||2^ and that the present data tend to prefer 
those parametrizations with large gluon content at small x (x ~ 0.1). 

The recent data on the jet and charm production in the two-photon process thus give us 
evidence and some quantitative information of the gluon content of the photon, but they are 
not yet accurate enough to determine directly the effective gluon distribution. It is therefore 
desirable to have a set of effective parton distributions in the photon with systematically 
different gluon distributions, so that we can learn more about the photon structure from 
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these experiments. 



In this paper we study all the available photon structure function data |TT|-p!^,p^-p6| at 
4 GeV^ < < 100 GeV^ in the leading order of perturbative QCD and, find a new set of 
the effective scale-dependent parton distributions in the photon, named WHITl to WHIT6, 
which are all consistent with the present data of the photon structure function and have 
systematically different gluon contents. We study carefully the charm quark contributions 
to the observed structure functions, which are evaluated by using the lowest order quark 
parton model matrix elements (7*7 cc and '~f*g — cc) and the massive Altarelli-Parisi 



(AP) equations We find that the photon structure function has a sensitivity to the gluon 
distribution at small x, but that a careful analysis is needed to determine experimentally 
the photon structure in this region. Predictions are also given for the total charm quark 
pair production cross section in the two-photon collision process at e~^e~ colliders. 

We note here that the next-to-leading order correction to the massless inhomogeneous 
AP equations has been known for a while P,|28|-pO|], and some phenomenological analyses 
1^-^ were performed at this level. Recently, the complete next-to-leading order correction 
has been obtained for the massive quark production process [^]. We work in the leading 
order of QCD, nevertheless, since errors in the experimental data as well as the theoretical 
uncertainties associated with the gluon contents of the photon are so large that the leading 
order approach is more suited to discuss them systematically. 

The paper is organized as follows. In Sec. 0, we review the AP equations that govern the 
effective parton distributions in the photon and introduce the notion of 'valence' and 'sea' 
components of the quark distributions. We discuss our parametrizations of the initial quark 
and gluon distributions and study charm contributions to the structure function. In Sec. |T|, 
we introduce all available photon structure function data, give a criterion to remove several 
low X experimental data, and then make a global fit of the initial light-quark distribution 
functions by using the leading order inhomogeneous AP equations. The fit is repeated 
by systematically changing the magnitude and the shape of the initial gluon distribution. 
The minimal of the fit as well as the distribution of the deviation of each data point 



from the best fit curve is examined carefully. Six effective parton distributions which have 
systematically different initial gluon distributions, WHITl to WHIT6, are then introduced 
and examined. In Sec. |rV|, we discuss effective heavy quark distributions in the photon by 
using the quark parton model (QPM) and the massive AP equations. In Sec. |V|, predictions 
are given for the total charmed particle production cross section in the two-photon process 



at e~^e~ colliders. Sec. ^ summarizes our findings. The details on the numerical methods 
that we use in order to solve the AP equation with and without charm quark mass effects 
are give in Appendix 0, and the parametrizations of our effective parton distributions in 
the photon are described in Appendix |B[ 

II. MODEL 

In this section we explain the framework adopted in this work to calculate the effective 
parton distribution functions in the photon and the photon structure function Q^). 

A. Inhomogeneous Altarelli Parisi equations 

In the and x region where the lightest n/-fiavor quarks can be produced, we have n/ 
quark distributions and Uf anti-quark distributions in addition to the gluon distribution in 
the photon. These are denoted by qi{x, Q"^), qi{x, Q"^) = 1 to nj), and g{x, Q"^) respectively. 
Apparently, the relation qi{x,Q^) = qi{x,Q'^) holds for each flavor. 

The evolution of these parton distributions is described by the inhomogeneous 



Altarelli-Parisi (AP) equations in the leading logarithmic approximation [^. For mass- 
less n/-fiavor case the AP equations can be written as follows: 



^^^^^ = ^elP,,{x) + \P,,{x) ® g.(a;, Q^) + P,,{x) ® g{x, Q')] , (2.1a) 

at /vr zvr L -' 

Uf 

2Pgg{x) (g) J2 Qii^^ Q^) + Pga^x^^i) ® 9{x, 



dg{x, g2) _ asiQ^) o d ^ „ , d „ ^ ^ r^^2^ ^2 lb) 



dt 27T L i=l 

where z = 1 to nj, t = logQ^/A^ with A being the QCD scale, and -Pj/s are the parton 
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splitting functions [^] whose explicit forms are given in Eq. (|A4]) of Appendix ^. The 
convolution integral is defined as a{x) b{x) = /J dy / y a(x / y)b{y) . 

As we show in the latter sections, the charm quark contribution to the photon structure 
function can be most conveniently calculated from the lowest order quark parton model 
processes (7*7 — > cc and Y9 ~^ cc) at present energies (Q^ ~ 100 GeV^), by excluding 
the charm quark distribution in the photon. Only at higher we introduce the effective 
charm quark distribution that evolves according to the massive n/ = 4 AP equations of 
Ref. [^. At very high Q^, we may neglect the charm quark mass and employ the massless 
inhomogeneous AP equations of Eq. (|2.1| ) with ^/=4. The matching of the quark parton 
model description with the effective heavy quark distribution in the massive AP equations 
is discussed in Sec. fV\ 

Hence in the analysis of the present structure function data that probe the photon 
structure up to ~ lOOGeV^, we introduce only three light-quark distributions (n/ = 3). 
Furthermore, in order to find a plausible parametrization of these quark distributions for the 
fit, we find it convenient to introduce the notion of 'valence' and 'sea' quark distributions. 
The 'valence' quarks are produced by the photon and the 'sea' quarks originate from the 
gluons. According to these definitions, the valence and sea quark distributions are written 
in terms of the original quark distributions g^'s: 



3 /p4\_ 2/ 2\ 

qseaix,Q') = 2j2 ) I ;\J q.{x,Q'), (2.2b) 

where (e^) = 2/9 and (e^) = 2/27 for nj = 3. Note that the singlet and non-singlet quark 
distributions, qs{x, Q"^) and qNs{x, Q"^), respectively, are related to our valence and sea quark 
distributions by 

qs{x,Q^) = 2j2q^{x,Q^), 

i=l 

= qy{x,Q'^) + qsea{x,Q'^) , (2.3a) 

q^six,Q') ^2j2[e^ - {e')]q,ix,Q') , 
1=1 



(2.3b) 



The photon structure function F2{x,Q'^) can be written in terms of qv{x,Q'^) and 
Q ) as 

F^{x,Q')^2xY.eM^,Q'), 

i=l 



X 



qv{x,Q ) + (e )qsea{x,Q' 



+ heavy quarks . 



(2.4) 



Heavy quark contributions will be discussed in subsection [11 C|. When we neglect small mass 



differences among the light three flavors, we can express the u, d and s distributions in terms 
of and qsea- 



u{x, Q^) = lqv{x, Q^) + ^qsea{x, Q^) , 

6 D 

d{x, Q^) = s{x, Q^) = —qvix, Q^) + -qsea{x, Q'^) . 



(2.5a) 
(2.5b) 



The AP equations of Eq. (|2.1| ) with = 3 can be rewritten in terms of the valence-quark, 
the sea-quark and the gluon distributions: 



dq„{x, 



a 



dt 2n'^'^''" 



{e')P,,ix) 



2tt 



-Pqq{x) ® qvix,Q'^ 



dqseaix,Q^) OsiQ"^ 



dt 27r 
dgix,Q') 



2 • 3Pgg{x) ® g{x, Q^) + Pgq{x) ® qsea{x, Q^) 
Pgq{x) ® [qv{x, Q"^) +qsea{x, g2)}+p,,(x,3) 



(2.6a) 
(2.6b) 
(2.6c) 



dt 27r 

It is now clearly seen that the valence-quarks are produced by the photon, the sea-quarks are 
produced by the gluon, while the gluon is produced by the valence-quarks, the sea-quarks 
and the gluon itself. Once a set of initial parton distributions at = Ql is given, we can 
predict the parton distributions at any Q'^{> Ql)hj solving the above equations numerically. 
The numerical methods which we use to solve these equations are explained in Appendix 0. 



B. Initial parton distributions 



To solve the AP equations of Eq. ( |2.6| ), we have to specify a set of initial parton distri- 
butions at = Qq- All the non-perturbative features of the photon structure are included 
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in these initial conditions. We use Qq = 4 GeV^ throughout our analysis in order that our 
perturbation approximation works well. 

As an initial valence-quark distribution, we take the following functional form: 

XQvix, Ql)/a = Aa;^"(l - xf^/B{B, + 1, + 1) , (2.7) 

where A^, By and are the free parameters which will be fitted to the experimental data, 
and B{a,j3) is the beta function that ensures the normalization, 

{xqy{x,Ql))/a= [ dxxqy{x,Ql)/a = Ay , (2.8) 
Jo 

for the energy fraction {xqy{x, Qq)) of the valence-quarks in the photon. 
As for the initial gluon distribution, we adopt the simple form 

xg{x, Ql)/a = Ag{Cg + 1)(1 - xf^ , (2.9) 

with two parameters, Ag and Cg. Again the normalization factor is chosen such that 

{xg{x,Ql))/a= f dxxg{x,Ql)/a = Ag (2.10) 

The present structure function data are not accurate enough to determine the gluon param- 
eters Ag and Cg. We therefore perform fit by the valence-quark parameters. Ay, By and 
Cy, by varying systematically the normalization {Ag) and the shape (Cg) of the initial gluon 
distribution. 

Before starting the fit to the data, we discuss plausible range of the gluon distribution 
parameters that we should explore. We obtain constraints on the ratio of the gluon energy 
fraction (Ag) to the valence-quark energy fraction (Ay) as follows 0. At sufficiently low 
momentum transfer scale (Q^ ~ 0.5 GeV^), only the long wavelength components of the 
photon are probed and the quark-antiquark pair produced from the photon undergoes non- 
perturbative soft QCD djTiamics that resembles the one which makes the quark-antiquark 
pair form the low-lying vector boson. The photon structure is then expected to have compo- 
nents similar to those of the vector bosons, in particular the p meson that couples strongly to 
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the photon. Although we do not know the structure of p, we expect the soft QCD dynamics 
to be insensitive to the total spin of the system and that it may be similar to the observed 



vr structure ||3^. If the photon had only this soft component, its deep structure should also 
be similar and we expect 

(x^(x,Qg))y ^ {xg{x,Ql))^ ^ ^ ^2 

In fact this ratio is common in the nucleon structure as well |35|j27[] and we can regard this 
ratio as an universal one from soft QCD dynamics. 

The photon, however, differs from the vector boson in that it is a source of a quark- 
pair with an arbitrary short wave-length. As the momentum transfer scale grows {Q^ ~ 
0.5 GeV^), one is more and more sensitive to these short wave-length components which 
dominate the photon structure at asymptotically high Q^. Although the transition from 
the regime where the vector meson-like component dominates to the regime where the short 
wave-length component dominates is gradual and it is governed by the non-perturbative 
dynamics, we may infer the effect of the latter component from its asymptotic behavior that 
can be calculated perturbatively. In particular, for the ratio of the gluon to the valence-quark 
energy fraction, we expect 



{xg{x,Q^)) 



{xq^{x,Q^)) 



1 (2,12) 
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for three light quark flavors [uf = 3). At the momentum transfer scale Ql = 4 GeV^, it is 
hence natural to expect the ratio to lie somewhere between the two extremes Eq. (|2.11|) and 



1 ^ {xg{x,Ql)) ^ ^ < , ^2 1 SI 

3 {xq^ix.QD) A ~ • ^ ■ ' 

We shall see in the next section that the valence-quark fraction is determined to be about 
unity by the experimental data of F2{x,Q'^). We will hence examine the parameter range 
1/3 ~ ~ 1 for the gluon energy fraction. 

Finally, we note that the sea-quark distribution is intimately related to the gluon distri- 
bution and that one cannot choose them independently. Although the sea-quark distribution 
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is in principle observable from the small x behavior of the photon structure function, we find 
that the present experimental determination of the small x structure of the photon suffers 
from an uncertainty associated with the unfolding technique adopted by most experiments 
p6| : this will be discussed in Sec. We therefore estimate the input sea-quark distribution 
by using the quark parton model cross section for the process '~f*g — > qq: 



xqsea{x,Ql) = 3^^ f dywl-,^] g{y,Ql), (2.14) 
Ztc Jax \y QqJ 

where a = 1 + irn^/Ql and 



w{z, r) = z 



P {-1 + 8;z(l - z) - Arz{l - z)} 
+ [z^ + (1 - zy + Arz{l - 3z) - 8r^z^] log 

1 fj 



(2.15) 



with P = yl — Arz/ (1 — z). The sea-quark mass m, which is taken to be common for the 
three light flavors, rriu = ma = = m, plays the role of the cut-off and we choose it to 
be 0.5 GeV. Here and throughout our analysis we adopt the leading order form of the QCD 
running coupling constant 

25. ..Q' 1,.. Q'^,^, 



with A4 = 0.4 GeV and = 5 GeV. Note that the effective number of quark flavors that 
governs the running of the coupling constant is chosen independently of the number n/ of 
massless quark flavors in the AP equations of Eq. ( |2.1|) . An accurate prescription for the 



choice of the effective number of flavors is found only in the next-to-leading order level 

We remark here that the above prescription leads naturally to a larger sea-quark input 
as the gluon input is enhanced. In particular, we find for the energy fraction ratio that 

holds almost independently of the input gluon parameters Ag and Cg in the region which 
we will discuss. The ratio increases with decreasing light-quark mass, and reaches 0.2 at 



m ~ 0.3 GeV. The ratio is about 0.3 in the parametrization of the vr structure |^4[, and its 
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perturbative asymptotic value is calculated to be 0.16 for nj = 3. Our sea-quark input is 
hence rather conservative for a given input gluon distribution. 



C. Charm contribution 

The charm quark cannot be incorporated into the massless AP equations in the region of 
moderate Q^, say < 100 GeV^, which has so far been probed by experiments. We should 
take into account the quark mass effect by using the massive-quark AP equation of M. Gliick, 
E. Hoffmann, and E. Reya and more accurately by incorporating the full next-to- leading 
order corrections We find by comparing with the results of the leading order massive- 
quark AP equations that the charm quark contribution to the photon structure function is 
well approximated by the sum of the contributions from the quark parton model processes 
7*7 cc and Yd cc at Q"^ < 100 GeV^. Beyond ~ lOOGeV^ the radiation of gluons off 
charm quarks is no longer negligible and, we should solve the massive-quark inhomogeneous 
AP equations. At large enough Q^, the charm quark mass effect to the Q^-evolution would 
become negligible and we can use the massless AP equations of Eq. ( p.l| ) with nj = 4. The 
matching of the distributions should then be made at appropriately large Q^. We will find 
that the charm-quark mass effects are not negligible even at ~ lOOGeV^. With the 
same criterion, the bottom quark contribution can be estimated by the lowest order process 
7*7 —>■ bb and '~f*g —>■ bb up to about ~ 1000 GeV^, above which we may introduce the 
effective b-quark distribution that follows the massive Uf = 5 AP equations. More accurate 
quantitative treatment [^] will become useful in the future when both the quark and gluon 



distributions are measured accurately from experiments. 

The charm-quark contributions to the photon structure function are thus calculated by 
the quark parton model at < lOOGeV^. The contribution of the direct process (7*7 cc) 
is given by 



a . ( ml' 



FUx,Q')Urect = 3-e>( X, ^ ) , (2.18) 
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where Cc = 2/3 is the charm-quark electric charge and the function w{x,r) is given in 
Eq. ( |2.15| ). In our numerical analysis, we take nic = 1.5 GeV. For the resolved process 
{l*g ^ cc) , we have 



FU^,Q%esoWe, = ^^^el / rfy^ -,;^ U(2/,Q^), (2.19) 



where a = 1 + Aml/Q'^, and the gluon distribution g{x, Q"^) is given by solving the massless 
= 3 AP equations of Eq. ( |2.6[ ) with the initial parton distributions of Eqs. (|2.7|), ( p.9| ) 
and 



The validity of our simple quark parton model calculation depends on how much the 
gluon emission by the charm quark distorts the effective charm quark distribution in the 
photon. The magnitude of this effect can be studied by using the massive AP equations for 
the charm quark and is presented in Sec. fTV] . 

III. FIT TO THE DATA 
A. Data 

In order to find good initial parton distributions at the energy scale Qo = 4 GeV^, we 
refer to all the available experimental data of the photon structure function at > Qq- In 
our analysis we use the data obtained by 8 groups at PETRA, PEP, TRISTAN and LEP 
e+e" colliders which are listed in Table |. 

We note here that not all the experimental data points are taken into account in our fit. 
First, we do not use the data at {Q"^) lower than 4.0 GeV^. Second, we accept only those 
data points at small x where the following inequality holds: 

^ {Q') .... 

^lower edge of the bm ^ ^q2^ _|_ ^\Yma.xy ■ \'^-^) 

Here is the experimental cut on the visible invariant mass of the final hadron system. 

Since those data points that violate the condition Eq.( p.lD are obtained at near the boundary 
of the experimental acceptance and since the sea-quark contribution to the structure function 
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can be rather singular at the low x region, they may suffer from large systematic uncertainties 
in the unfolding procedure ||36|. As a result of the above two requirements, 47 data points 
are retained in our fitting, which are all listed in Table |. 



B. Fit 

By fitting our theoretical predictions for the photon structure function to these exper- 
imental data, we tune the parameters of the initial valence-quark distribution A^, and 
C^. We repeat the fit by varying the initial gluon distribution parameters Ag and Cg sys- 
tematically, while keeping the strong coupling constant and the charm quark mass fixed at 
A4 = 0.4GeV and rric = l.SGeV, respectively. In particular, we examine the case with Ag 
= 0.5, 1 and 1.5 systematically by changing Cg and find little sensitivity of the structure 
function data to the shape parameter Cg. The fit results for arbitrarily chosen 12 cases 
are summarized in Table 0. Although there is a tendency in the data that prefers small 
Ag (small gluonic energy fraction) and small Cg (hard gluon distribution), it is caused by 
a few data points at small x with relatively large deviations from the best fit curve, as we 
will discuss below. We choose three representative Cg values {Cg=3, 9, 15) for each of the 
normalizations Ag=0.5 and 1 that are consistent with the ansatz Eq. (|2.13|) . These gluon 
inputs are named as WHITl to WHIT6, respectively, as shown in Table |nT[ 

Fig illustrates the matching of the data and the theoretical curves, and Fig. |^ shows 
the distribution of the deviation of each data point from the best fit value, {F2{x)fit — 
i^2(x)data)/t^(i^2(a;)data). As cau be seen from Fig. | and Fig. |, all of WHITl to WHIT6 
gluon distributions give similar quality of fits to the photon structure function data. The mild 
dependence of the best fit value on the choice of the initial gluon distribution parameters 
as shown in Table ^ is a consequence of a few data points at the lowest x bin that satisfy the 
criterion Eq. ( p. 11 ): see Fig. 0. In Fig. ^ deviations of those data points that are removed 
from the fit by the criterion Eq. ( |3.1| ) are also indicated by dashed lines. The large Ag large 
Cg gluon distributions lead to a significant rise in the structure function at small and 
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a naive integration of the structure function in the given x-bin tends to be very sensitive 
to the lower edge of the lowest x-bin as is seen in Fig. |1]. After imposing the selection 
criterion Eq. ( |3 . 1| ) , such sensitivity to the very low x behavior of the structure function 
is almost completely lost, as can be seen from the deviations of the data points that are 
connected by the solid lines in Fig. 0. In view of the relatively large theoretical uncertainty 
in simulating hadronic events at small x, we conclude that the present experimental data 
on the photon structure function have poor sensitivity to the gluonic content of the photon. 
The normalization of the valence-quark distribution Ay is found to be roughly 1, regardless 
of the difference in the sea-quark contribution that depends strongly on our gluon inputs. 

We find from Table that the best fit values of the initial valence-quark parameters are 
almost the same for different C^'s for a common Ag. Hence we introduce a 'standard' set of 
the valence-quark parameters for each Ag] i.e. 

Ay = 0.94 , By = 0.50 , Cy = 0.25 for Ag = 0.5 , (3.2a) 
Ay = 0.89 , By = 0.70 , Cy = 0.45 for Ag = 1.0 . (3.2b) 

We calculate the values for various C^'s by fixing the normalization Ag and the associated 
valence-quark inputs as above, and the result is summarized in Table The Cg dependence 
of the resulting is also presented in Fig. The minimal values as obtained in Table |T| 
by tuning the valence quark parameters are also shown by large symbols. The tuning of the 
valence-quark parameters do not improve the fit much. We therefore adopt the common 
valence-quark input of Eq. (|]2|) for the sets WHITl to WHIT3, while that of Eq.(|3]2bD 
for the sets WHIT4 to WHIT6, and hereafter we call the parton distributions with these 
inputs as WHITl to WHIT6. The slightly small valence-quark contributions of Eq. ( |3.2bD 
at small x partially compensates for the larger sea-quark contributions associated with the 
large gluon inputs of WHIT4 to WHIT6. 
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C. Gluon distribution 



As we described above, we present a set of six effective parton distributions in the photon 
with systematically different gluon contents. We show in Fig. ^ all the gluon distributions 
(WHITl to WHIT6) at three momentum transfer scales, (5^=4, 20 and 100 GeV^. The 
area under the curves at Q^= 4GeV^ is given by the normalization Ag = 0.5 for WHITl to 
WHITS, and = 1.0 for WHIT4 to WHIT6. The shape of the distribution becomes softer 
as we move from WHITl to 3, and from WHIT4 to 6 in each set. The huge difference in the 
initial gluon distributions tends to diminish at higher Q^, as expected from the asymptotic 
behavior of the solution of the inhomogeneous AP equations. 

Also shown in Fig. ^ for comparison are the effective gluon distributions of GRV [^, 
DG and LACl ||10|] at the three momentum transfer scales. We note that in the x 
and range as shown in Fig. ^, our WHITl effective gluon distribution is similar to the 
gluon distribution of GRV p||, while WHIT6 gluon distribution behaves similarly to that of 



LACl |T0|. It should be noted, however, that LACl parametrization for the effective gluon 
distribution is more singular at x — > than that of WHIT6, which results in a very large 
energy fraction {xg{x, Ql))/a = 2.37 at Qo = 4 GeV^, in conflict with our ansatz Eq. (p.l3|) 



with ~ 1. Accurate measurements of the photon structure functions at very small x 
as well as the high energy behavior of the charm and mini-jet production cross section at 
77 and 7p collisions will be able to distinguish the different small x behavior of these two 
effective gluon distributions. 

IV. EFFECTIVE HEAVY-QUARK DISTRIBUTIONS IN THE PHOTON 

In this section, we compare the result of the quark parton model (QPM) calculation of 
the effective charm quark distribution with that of the massive inhomogeneous Altarelli- 
Parisi(AP) equations |2^. We expect that the QPM approach is appropriate at low mo- 



mentum transfer scale where the charm quark mass effect is significant. At high momentum 
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transfer Q"^ /ml S> 1, gluon emission from charm quarks is no longer negligible and we need 
the massive-quark AP equations to sum up the leading effects. On the other hand, the ap- 
proach of neglects the charm quark mass effects in gluon emission from charm quarks, 
and hence it may overestimate the gluon emission effects at low . We therefore use the 
QPM prescription for the effective charm quark distribution at < 100 GeV^ and switch 
to the solution of the massive quark AP equations at > 100 GeV^. 

The QPM charm quark distribution consists of the valence part and the sea part and is 
defined as 



.QPM 



,QPM/ 



(4.1) 



where 



.QPM, 



2xe1 



2xel 



-^2^c('^5 Q ) [direct ) 
P'2,ci.'^yQ )|resolved) 



(4.2a) 
(4.2b) 



with F^^{x, 



direct, 



resolved as givcu in Eq.( p.l8| ) and Eq.( p.l9| ), respectively. 



As is clear from these definitions, the effective charm quark distribution calculated by 
the QPM reproduces the photon structure function well, but we should expect large process- 
dependent threshold corrections when it is used for other processes with an equivalent real 
charm quark in the photon. 

At high Q^, we expect that emission of collinear gluons from charm quarks becomes 
significant and we solve the massive-quark inhomogeneous AP equations 



dqi{x,Q'^] 
dt 

dg{x,Q'^) _ a^iQ'^ 

dt ~ 
dc{x, Q"^) 



[Pgqix) ® qi{x, Q2) + p^^^x) ® g{x, Q^) 



2n 



dt 



a 
2^ 



Pg,{x)®2 ( Y.^i{x, Q')+c(x, Q2) j +Pgg{x, A)®g{x, Q^) 
Pqg{x) ®c(x, Q^)+Pcg{x, Q2) ®g{x, Q^] 



„ (ic'Pc'yi^XjQ ) 



27r 



(4.3a) 
(4.3b) 

(4.3c) 



where i = u, d, s, and and Peg are the photon and gluon to massive-quark splitting 
functions, respectively, as defined by ^7 
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/ at X 



with a = 1 + Aml/Q'^. The function i(7(x,r) is defined in Eq. ( p.l5| ). 



As in the case of hght quarks, it is convenient to divide the charm-quark distribution 
into the valence and the sea parts: 

c(x, Q2) = c„(x, Q^) + Cseaix, Q^) . (4.5) 

The valence-charm-quark comes from the photon and the sea-charm-quark comes from the 
gluon. Eq. (|4.3c| ) can then be split into the following two equations for and Csea'- 



eiPc^ix, Q') + ^^P,,(x) ® c,(x, Q') , (4.6a) 

Pm{x) ® Csea{x, Q') + Pcg{x, Q^) ® 9{x, Q^)] , (4.6b) 



dt 27r ^ ' ' 27r 

c?Csea(a;, Q^) «s(Q^ 
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The massive-quark splitting functions Pcg{x, Q"^) and ^£7(2;, Q^) are singular at the charm 
threshold x = 1/a. Due to this singularity, we find that much CPU time is needed in order 
to get an accurate numerical solution when one solves Eqs. ([4.3a]) , (|4.3b|) , ( |4.6a| ) and ( |4.6b|) 



directly. The numerical problem associated with the use of the massive-quark splitting 
function of is severer in the photon structure than in the proton structure because of 
the presence of the leading inhomogeneous term. We can avoid, however, appearance of 
the singular massive splitting functions, and obtain a set of equations which contain only 
smooth functions by dividing further the valence- and sea- charm-quark distributions into 
the QPM part and the remnant. The QPM part of the valence- charm- quark distribution is 
defined by Eq.( 4.2a|) and that of the sea-charm-quark distribution is defined by Eq.( 4.2b| ). 



The remnants are defined through the following equations: 



c„(x, Q2) ^ cQ™(x, Q2) + 5c,(a;, Q^) , (4.7a) 

sea{x, Q') = cZ^'ix, Q') + Sc^U^, Q') . (4.7b) 



Using Eqs.(|4.2|), (|4.4|), and ([4.6|), one can derive the following equations that govern the 



deviations 5c„ and Scgea'- 
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dt"^"^"'^ , 27r 



5c,ea(x, = ^:^^P,,(x) ® cQP^(x, Q^) + 5c,,,(x, Q') , (4.8b) 



within the similar approximation that Gliick, Hoffmann and Reya made ||2^. Here we take 
the boundary conditions 

Sc,,sea{^,Ql)=0. (4.9) 

Note that there now appear no massive sphtting functions in these equations. Intuitively, this 
is because, the deviation of the effective charm quark distribution from the QPM prediction 
is caused by the gluon emission from the charm quarks, and the emission is approximated by 



the massless quark splitting function Pqq{x) in the scheme of Ref. ||2^. A brief explanation 
of the numerical method which we employ to solve these equations is given in Appendix 0. 

The resulting effective charm-quark distributions (multiplied by x) in the photon are 
shown in Fig. |^ for (5^=4, 20 and 100 GeV^, where the boundary conditions ( |4.9| ) are set at 
Ql =4 GeV^. Also shown in the figure for comparison are the valence-up- quark distributions 
of WHIT 1-3 and WHIT4-6. Predictions of the QPM are shown by solid lines while those of 
the massive-quark inhomogeneous AP equations are shown by dash-dotted lines. We find 
that the differences between the predictions of the two approaches are negligibly small for 
the sea-charm-quark distribution in WHITl to WHIT6. From Fig. |^, we find that the 
QPM prediction for the valence-charm-quark distribution differs by up to about 20% at 

= lOOGeV^. The shape of the prediction in the massive inhomogeneous AP equations is 
softer than the QPM one as expected. Since the massless splitting function Pgg is used for 
describing the gluon emission from the charm quark, the deviation from the QPM predictions 
in the lower region may be an overestimate. 

From the above discussion, we conclude that the QPM calculation of the effective charm 
quark distribution is appropriate in the region ^ lOOGeV^, given the present experi- 
mental accuracy and theoretical uncertainties in the charm quark mass and higher order 
corrections. At sufficiently high Q^, the massless 4-flavor AP equations should become a 



good approximation. However, we find that the charm quark threshold effect is still signifi- 
cant near x ~ 1 even at ~ 100 GeV^. We therefore use the massive inhomogeneous AP 
equations between = lOOGeV^ and = 2500GeV^ up to where we give parametrizations 
of the parton distributions. We require the continuity of all effective parton distributions at 
Q'^ =100 GeV^: at 4 GeV^ < < 100 GeV^ we use the massless n/ = 3 inhomogeneous 
AP equations with the QPM approximation to the charm quark distribution, and at 100 
GeV^ <Q^< 2500 GeV^ we use the massive inhomogeneous AP equations with the massive 
charm quark. The boundary conditions (|4.9| ) are hence set at Qq = 100 GeV^. 

Finally, as is mentioned in Sec. [II C| , the effective bottom quark distribution in the photon 
may be approximated by the QPM calculations (7*7 — > bb and '~f*g bb) all the way up 
to ~ 2500 GeV^, and hence we do not give parametrizations of the effective b-quark 
distribution in the photon. 



V. CHARM PRODUCTION CROSS SECTION 



In this section we study the charm-quark production cross section via the two-photon 
processes by using our new effective parton distribution functions in the photon. The charm 
quark production cross section is expected to be much more sensitive on the gluonic content 



of the photon than the photon structure function [23 



A. Equivalent real photon approximation 

To calculate the charm-quark production cross section in the two-photon processes, we 
employ the equivalent real photon approximation (EPA) to the nearly on-shell virtual pho- 
tons 1,0. 

In the EPA the charm-quark pair production cross section for the process e'^e~ — *• 
e~^e~ccX is approximated by the convolution of the effective real photon fluxes in the 
beam and the cross section of the subprocess 77 — > ccX: 
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a{e~^e — > e"*"e ccX) = J dxidx2D^/e{xi, Qf)D^/e{x2, Ql) d-{'ji'j2 ccX) , (5.1) 

where (t(7i72 —* ccX) is the subprocess cross section at = sxiX2, and D^f/e{xi, Qf) {i = 
1, 2) is the equivalent real photon distribution in the electron {i = 1) or in the positron(i = 2). 
The cm. energy of the colliding e~e"*" is -y/i, and that of the colliding two photons is W. 



The improved form of the photon distribution is written as |^ 

a l + (l-x? 



D-y/e{x,Q 



log^-1 



^min 



+ ^x, (5.2) 
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where tmin = mlx"^ / (1 — x) is the kinematical limit of the magnitude of the lepton momentum 
transfer t and the second term in the right-hand-side denotes the contribution of the electron 
helicity flip amphtudes ISM]. The scale should be set by the dynamical condition that 



the subprocess "cross section" for a virtual process is damped at |t| > 

We choose the scales Qf as follows, when the photon (74) couples directly to the cham 
quark, we set 

Ql = min[m^ -hp4,tmax,tcut], (5.3) 

where pt is the transverse momentum of the charm quark in the 77 cm. system and tmax = 
s{l—Xi) is the kinematical maximum of the momentum transfer, and tcut denotes the possible 
experimental cut on the magnitude of the momentum transfer. The scale + can be 
interpreted as the virtuality of the internal charm quark line; if the virtuality of the photon 
is larger than that of the internal charm quark, the production of charm quark is strongly 
suppressed. If the photon 7j resolves into light partons that contribute to the cc production 
subpro cesses, we set 

Q2 = min[2GeV2, t,^,], (5.4) 

since the gluon content of the photon should be suppressed if its virtuality is much larger 
than the hadronic scale. Although in principle the quark and gluon contents of a sufficiently 
virtual photon is calculable in perturbative QCD, we neglect contributions from virtual 
photons with >2 GeV^. The uncertainty associated with the choice of the cut-off scale 
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may be estimated by changing the scale by a factor of two, and it is at a few % level at 
TRISTAN/LEP energies, while it can be reduced significantly by reducing the tcut value by 
an anti-tagging. 

To check the validity of our approximation Eq. ( |5.1| ) with the scale choice Eq. (|5.3| ), we 
compare the EPA prediction with the exact cross section for the direct process e+e^ 
e~^e~cc in which both photons couple directly to the charm quarks. In Fig. ^ we compare 
differential cross section of charm quark pair production evaluated exactly (solid lines with 
error bars) with that of the EPA prediction (dotted lines) against the cosine of the scattering 
angle 6 of the subprocess 7*7* —>■ cc in the cc rest frame. The exact cross section has been 



calculated by HELAS [p9| with the help of a Monte Carlo integration package BASES [40 



Near the charm quark pair threshold (Fig. |^a) we find almost exact agreements between 
the EPA and the exact results. At far above the threshold energy (Fig. ^jb), the subprocess 
cross section has peaks at cos 6 = ±1 where our EPA prediction underestimates the cross 
section by ~ 30%. Fig. ^ shows the charm quark pair invariant mass ("\/S) distribution 
after integrating over the scattering angle 6. Though there is a tendency that our EPA gives 
slightly smaller cross section, we find that the EPA calculation of the total cross section 
agrees with the exact cross section at 1 % level. 



B. Charmed particle production cross section 

The inclusive 77 — ^ ccX cross section in Eq. (|5.1|) , including the resolved processes, is 
described in the leading order as 

d(T(77 ccX) = (io-(77 cc)\s=w^ 

+ 2 J dx g{x,Q'^)da{jg cc)\s=w^x 

+ J dxdy g{x, Q'^)g{y, Q"^) da{gg cc)\s=w2xy 

+ 2 dxdy (li{x,Q'^)qi{y,Q'^) da{qq ^ cc)\s=w^r,y , (5.5) 
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where W is the invariant mass of the two photons and is that of the charm quark 
pair. The effective gluon and quark distribution functions in the photon are denoted by 
g{x,Q^), QilXjQ"^), for which we use the parametrizations of WHITl to WHIT6 parton 
distributions. The scale of the parton distribution functions and the QCD couphng 
as{Q^) has been chosen as follows: Q"^ = ml + for •yg — >• cc and gg cc, and = s for 
qq — > cc. 

We show in Fig. |^ the total charmed particle production cross sections. The six curves 
are obtained by setting = l.SGeV, a = 1/137 and A4 = 0.4GeV, and by imposing 
the open charm production cut W > 2m£, = 3.74 GeV. The vertical bars attached to the 
WHITl, WHIT4 and WHIT6 curves indicate typical uncertainties in the cross section on 
the charm quark mass between 1.3 GeV and 1.7 GeV. It is notable that uncertainty due 
to the charm quark mass is not reduced at high energies, because the total cross section 
is always dominated by the charm-quark pair production near the threshold. Higher order 
QCD corrections also modify the cross sections. 

As can be seen from Fig. |, the total charm quark production cross sections at high 
energy e"'"e~ collider experiments are sensitive to the gluon distribution in the photon, which 
cannot be measured accurately at the current photon structure function experiments as we 



describe in Sec. |T|. At TRISTAN energies, the predictions using WHITl to WHIT3 and 
those using WHIT4 to WHIT6 are almost the same within each group. At LEP and LEP200 
energies, the total charmed particle pair production cross section of the two photon process 
is more sensitive to the small x behavior of the gluon distribution, and WHITl to WHIT6 
predictions can be distinguished provided that the uncertainties due to the charm quark 
mass and higher order corrections are reduced. 

Fig. H shows the total cross sections for the process e+e^ e^e^ccX, together with 
the individual contribution of the direct photon process (the first term of Eq. (p.5|)), the 
once-resolved photon processes (the second term of Eq. ( p.5| )), and the twice-resolved photon 
processes (the third and the last terms of Eq. ( |5.5| )). From Fig. ^ we find that the resolved 
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processes, which are governed mainly by the gluon contents of the photon, grows much more 
rapidly for WHIT6 than for WHITl. The contribution from the once-resolved processes 
overcomes that of the direct process even at TRISTAN energies for WHIT6, while they be- 
come comparable only at around LEP200 energy for WHITl. Hence the energy dependence 
of the once-resolved photon contribution to the charm quark pair production cross section 
will be useful to determine the gluon distribution. 

VI. CONCLUSIONS 

We have studied all the available photon structure function data [p!T| - [T7| , p^ -p^ at 4 
GeV^ < < 100 GeV^ in the leading order of perturbative QCD, and have found a new 
set of the effective scale-dependent parton distributions in the photon, WHITl to WHIT6, 
which are all consistent with the present data (Fig. |1]). The six parton distributions have sys- 
tematically different gluon contents (Fig. ^) and their parametrizations are given in Table ^ 
to |X[ We have studied carefully the charm quark contributions to the observed structure 
functions, which are evaluated by using the lowest order quark parton model matrix ele- 
ments (7*7 — i> cc and 7*5' — cc), and by using the massive inhomogeneous AP equations. 
We have found that the photon structure function has poor sensitivity to the gluon distri- 
bution, except at very small x. In order to probe the gluon content of the photon from the 
photon structure function at small x, a careful analysis of experimental data is needed. 

Predictions have also been given for the total charm quark pair production cross section in 
the two-photon collision process including the resolved processes at e^e~ colliders. At PEP 
and PETRA energies, the difference in the predictions of WHITl to WHIT6 distributions 
is badly observable, while at TRISTAN energies WHIT4 to WHIT6 distributions predict 
significantly higher cross section than WHITl to WHITS (Fig. For all of our parton 
distributions, WHITl to WHIT6, the contribution of the once-resolved process exceeds that 
of the direct process at energies above around 200 GeV (Fig. ^. 
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APPENDIX A: THREE METHODS TO SOLVE THE AP EQUATIONS 

In this appendix we introduce three numerical methods for solving the massless inhomo- 
geneous AP equations of Eq. ( |2.1| ) and the massive inhomogeneous AP equations of Eq. ( [4.3[ ). 
We first introduce the Mellin transformation technique in Sec. |A 1| , and in Sec. [A2| we de- 



scribe the recursion method that we actually use in performing the fit to the experimental 
data. We use the Runge-Kutta method to solve the massive inhomogeneous AP equations. 



which is discussed in Sec. |S~3 



The first two methods, the Mellin transformation and the recursion method, work only 
for the massless AP equations, while the Runge-Kutta method can be applicable to both 
the massive AP equations and the massless ones. We use all these three methods for cross 
checking of our numerical results for the massless AP equations. 

1. Moment method for the massless inhomogeneous AP equations 

The Mellin transformation jin) of ][x) is defined as 

~f{n)^ f'dxx'^-'fix) , (Al) 
Jo 
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where n is a complex number. This transformation solves the convolution integrals in the 
inhomogeneous AP equations as simple products of the Mellin transforms of the splitting 
function and the parton distribution function: 

r dxx^-^Pix) ® qix) = P(n)q(n) . (A2) 
Jo 

The inhomogeneous AP equations with Uf massless quarks 

= ^e^P,,ix) + [P,,ix) q,{x, Q') + P,,{x) ® g{x, Q^) 

dg{x,Q') a,(Q2) 



dt 27r 
with the splitting functions 



Pgq{x)^2Y,qi{x,Q^) + Pgg{x,nf)^g{x,Q^ 



Pq-y{z 



1 + '^' 3^,^ 
0^ + 2^^'-^^^. 



4 1 + 

3 z 

^ V + (1 - 



(l-z). 



11 



12 18/ ^ 



(A3a) 
(A3b) 



(A4a) 

(A4b) 
(A4c) 

(A4d) 
(A4e) 



are transformed to the following form in the n-space: 

'Pggin)qi{n, Q^) + Pgginyg{n, Q^) 



with 



dqi{n, Q^) _ ^ 2p / N , o^siQ^) 
dt -27r''' "^^""^^ 27r 



dg{n,Q'') as{Q^) 



dt 



2n 



nf 



2 E PMq,{n, Q') + P,,{n, nf)g{n, Q') 



(A5a) 
(A5b) 



Pqqin 



Pgqin 



Pgg{n 



Pggin 



Pqj{n 



3 _ 2n + 1 

4 ~ 2n(n + 1) 
An^ + n + 2 

3 n{n^ - 1) ' 
n'^ + n + 2 
2n{n + l){n + 2) ' 

6 



- 7 - 



11 ^/_|_ 2 — n 



12 18 n{n-l) (n+l)(n + 2) 



- 7 - ^'(n) 



6P,,{n) 



(A6a) 

(A6b) 

(A6c) 

(A6d) 
(A6e) 
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Here = d\ogr{x)/dx is the di-gamma function and 7 = 0.577- ■ ■ is Euler constant. 

Eq. ( |A5| ) can be readily solved analytically by diagonalizing with respect to the parton 
flavors, and we flnd qi{n, Q"^) and g{n, Q^) at arbitrary > Qq. 

The x-space solution, qi{x,Q'^), can then be obtained by performing the inverse Mellin 
transformation 

g,(a;,Q2) = / rfnx-"g,(n, Q^), (A7) 
2m Jc 

numerically. The complex integration path C must be in the right half plane of all singu- 
larities of the integrands. We choose the path 



C = C+ - C± : n = 2.5 + exp(±37rz/4)M, {u goes from to cx)) . (A8) 

The angle of this path, 37r/4, makes the integrals converge. For each value of x in Eq. ( [A 71 ), 
\x~^\ = I exp(— log(a:;))| ~ exp(log(a;)M/v^) can be neglected at sufficiently large u e.g. 
u ^ — 25/log(a;), for slowly varying n-space functions, qi{n,Q^). 

Although the Mellin transformation method is compact and fast, it is not generally useful, 
because we have to restrict our input distributions qi{x,Ql), g{x,Ql) to those functions 
whose Mellin transformations can be analytically obtained and hence we can not use it to 
solve for our sea-quark distributions (see Eq.( |2.14 )), and once the splitting functions Pij{n) 



have a certain dependence on Q^, which is the case in the massive AP equations, we can no 
longer solve Eq. (^^) analytically in general. 



2. Recursive method for the massless AP equations 

The second method to solve the massless AP equations is based on the power expansion 
of the solution by the double logarithmic energy scale parameter s: 

s ^ log(t Ao) = log[ log(QVA')/ log(Q^/A') ] . (A9) 

To improve the convergence of the power expansion and the behavior at small x, we introduce 
the rescaled parton distribution gj, defined as 
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qi{x,s) = x'^e ^qi{x,Q'^] 
g{x,s) = x'^e~'g{x,Q'^) 



(AlOa) 
(AlOb) 



The massless inhomogeneous AP equations for Uf quark flavors, Eq. (|2.1| ), can then be 
rewritten as follows: 



dqi{x,s) 

ds 
dgjx, s) 

ds 



a 



io^e-Pg^(x)+toe' 



2n 



Pqq{x)®qi{x, S)+Pqg{x)®g{x, s) 



toe' 



2tt 



i=l 



-qi[x, s) 
g[x,s) , 



where we define the rescaled splitting functions, 



Piji^X^ X Piji^X^ , 



(Alia) 
(Allb) 

(A12) 



which is regular at x 0. Since the QCD running coupling ^^(Q^) scales as exp(— s) ~ 
1/ logQ^ in the leading log approximation, the factor e^as{Q'^) is regarded to be a constant: 

e = toe-^ = ' , (A13) 
" 27r 33 - 2ni. ' ^ ' 

where is the number of light- quark flavors which governs the running of the QCD coupling 

a,. We adopt n^. = 4 for 4.0 GeV^ < Q"^ < 100 GeV^, and nj. = 5 for Q"^ > 100 GeV^, 

according to Eq. (|2.16| ). 

Integrating Eq. (|A11| ), one finds 

qi{x,s) = qi{x,0) + sto—eiPi^{x) + eJ^Pij ^ ds'qj{x,s')- ds'qi{x,s'), (A14) 

Ztt j Jo Jo 

where qi{x, 0) is nothing but the rescaled input parton distribution at = Qq. Expanding 

the parton distributions by powers of s, i.e., 



qi{x,s} 



e=o 



X) 



one can reduce Eq. ( |A14| ) order by order of s: 

ql°\x) = qi{x,0) , 



g\x) = eY.P^A^)®qr'\x)-c^f-'\x) 



t,^ejP,^{x) + eY,h,{x)®qf\x)-(i\x) 
27r ^ 



5(^-1), 



(for £ > 2) . 



(A15) 

(Al6a) 
(A16b) 

(Al6c) 
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The above equations give the input parton distribution qi{x, 0) as the zero-th order approxi- 
mation for qi{x, s). The first correction which is hnear to s is the sum of the inhomogeneous 
term and the terms which is driven by the zero-th approximation (|A16b|) . The higher order 



corrections are determined recursively by Eq.( [A16cD . Summing all the contributions as in 



Eq. ( [A15| ) up to an appropriate order, one obtains the solution qi{x, s) with a given accuracy. 

This method is useful for arbitrary input parton distributions. This is an advantage 
of this method as compared to the moment method in Sec. |A 1|, since our input sea-quark 
distribution is calculated by the convolution integral in Eq. ( p. 141) and hence does not have 



analytic Mellin transform. Furthermore, we find that the recursive method needs much less 
CPU time than the more general Runge-Kutta method of Sec. |A 3| . We therefore use this 
scheme in the actual fitting of the experimental data. 

When one solves the homogeneous AP equations, similar method can be used, where 
the scale parameter s is expressed often by the QCD coupling constant as(Q^) as s = 
\og[as{Ql) / as{Q'^)]. In this definition, one can naturally incorporate the effect of the change 
of the effective number of the quark flavors np at the quark mass thresholds by simply 
rescaling the s variable. However, due to the presence of the inhomogeneous term, we 
cannot absorb all scales into the s variable. 

Even though the recursive method is powerful for solving the massless inhomogeneous 
AP equations, it cannot be used to solve the massive AP equations. This is because the 
signiflcant threshold effect in the massive quark distribution and the singular mass effect in 
the massive splitting functions do not allow the power expansion like Eq. ( |A15D 



3. Runge Kutta method 

The Runge-Kutta method is the general method to solve the differential equations and 
it requires a relatively large CPU time and sufficiently large number of data points for 
precise calculation. We use this method to integrate the inhomogeneous AP equations with 
a massive charm quark: 
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dqi{x,Q'^) 
dt 

dg{x,Q'^) 
dt 

dc{x, Q"^) 
Jt 



2^ i ) ^ 27]- qq\-^ ) Ht\-^i ^ J < qg\ 



[Pqgix) (g) + p^^^x) (g) g^x, Q^) 



27r 



ii=u,d,s 



^elP,,{x,Q^) + ^ 



Pgg{x)(gc{x, Q^)+Pcg{x, (g^g^x, Q^) 



(A17a) 
(Al7b) 

(A17c) 



After introducing the QPM component of the charm quark distribution by Eq. ( |4.2| ), the 
deviations from the QPM predictions as defined in Eq. ( |4.?| ) satisfy the following equations: 



cf"^{x,Q^) + 5c,{x,Q' 



dt 



(Al8a) 
(A18b) 

It is the above two equations with the boundary condition Eq. (|4.9| ) that we solve by using 

the Runge-Kutta method. It is difficult to apply the previous two methods, the Mellin 

transformation method or the recursion method, to solve the above equations. For the 

Mellin transformation method, we can not analytically Mellin transform c^fi^. For the 

recursive method of Sec. [A 2| , the expansion of w{z, r) in powers of s fails at the threshold 
^ = l/a= l/(l + 4m2/Q2). 

To solve the massive AP equations of Eqs. ( [A17a| ), ( [A17b| ), ( [Al^ ) we use the fourth 
order Runge-Kutta method with adaptive step size control for calculating the evolution of 
t. By discretizing the x variable as xj {j = 0, 1 ■ ■ ■ , A^), we solve the AP equations as a 
set of 5(A^ + 1) ordinary differential equations, while performing the convolution integral 
P{x) q{x) by the Simpson integration. We choose 



tanh (sinh(-2 + 0.04j)) /2 + 1/2, (j = 0, ■ ■ ■ , 110) , 



Xiii = 1 



(Al9a) 
(Al9b) 



so that the data points are dense near the both ends a; ~ and a: ~ 1. We find that A^ = 111 
is sufficient to archive an accuracy of 1%. 
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APPENDIX B: WHITl WHIT6 PARAMETRIZATIONS 



For practical use of the WHIT parton distributions which are given by the standard 
valence-quark input parameters of Eq. (|3.2| ) and the various gluon input parameters as 
listed in Table [TTTI , we give their convenient parametrizations (or simple prescriptions to 
calculate them). The parametrizations are given in three different regions, 4GeV^ < 
g2 < lOOGeV^, lOOGeV^ < < 2500GeV2, and < 4GeV2. The FORTRAN code 
for the distributions that we use for generating curves in this paper is obtainable from 
kaoru@kekvax.kek.jp or from kaoru@jpnkekvx.bitnet 



1. Parametrizations in the region 4GeV^ < < lOOGeV^ 

In this region, we parametrize those solutions of the massless 3-flavor AP equations 
which are described in Sec. |111 B|. The heavy quark distributions are calculated by the quark 



parton model as is mentioned in Sec. |IV|. Therefore, it is sufficient to give parametrizations 
for the valence-quark, the gluon and the sea-quark distributions. 

The valence-quark distribution is parametrized in the following functional form: 



xq^{x, s)/a = {Aq + Aix + AgX^) x^{l - x)'^ , (Bl) 

where s is defined in Eq. (^) with Ql = AGeV^ and A = A4 = 400MeV, and A^'s, B and 
C are polynomials of s of at most 4-th order. The coefficients of these polynomials can be 
found in Table fV| to Note that the valence-quark distributions are common for WHITl 



to WHITS and for WHIT4 to WHIT6. 

For the gluon distribution, we take the following form: 

xg{x, s)/a = Ax^{l - xf + {A'q + A[x) x^' {1 - xf . (B2) 

As in the case of the valence-quark, A's, i?'s, and C's are polynomials of s of at most 4-th 
order. The coefficients of these polynomials are also found in Table to |X[ Note that the 



second term in Eq. Q is common for WHITl to WHITS and for WHIT4 to WHIT6. The 
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reason is that it approximately represents the contribution of the gluons emitted from the 
valence-quark, and that the valence-quark distribution is common for WHITl to WHIT3 
and for WH1T4 to WH1T6. 

Since the input valence-quark distributions are much harder than the input gluon distri- 
butions for the WHIT parton distributions as seen in Eq. ( |3.2| ) and Table the contribution 
of the common part of the gluon distributions is expected to be significant at x ~ 1. In 
fact, we can estimate the behavior of the common part from Eq. ( p.6c|) . In the one gluon 



emission approximation, we get the following relation for the common gluon distribution 
which originates from the valence-quark: 

g{x) ~ Pgq{x) (g) q^{x) . (B3) 

Assuming the behavior of the valence-quark distribution at x ~ 1 as ~ (1 — x)*"", 

Eq.(g) leads to 

g{x) ~ (1 - x)^"+^ as X ^1. (B4) 

This common part of the gluon distributions is dominant at x ~ 1 in the WHIT parton 
distributions, because < 1 and Cg > 3 for all the input parton distributions. One 
may find that the above crude estimation of the behavior of the common part of the gluon 
distributions at x ~ 1 works rather well even in quantitative sense by looking at the relevant 
entries of Table to II^. 

As for the sea-quark distribution, we parametrize it by the following functional form: 



XQs 



,(x,s)/a = Ax(^o+^^")(l-a;)^. (B5) 



The coefficients that describe the s dependence of A, i?'s and C are also listed in Table [TVI 

torn 



The effective charm quark distribution is calculated by the quark parton model by using 
Eqs. ( [4.1[ ) and ( [4. 21 ). In the calculation of the sea-charm-quark distribution in Eq. ( |4.2b| ), the 



above parametrization of the gluon distribution for each WHIT parton distribution is used. 
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An efficient integration routine for the convolution in the calculation of the sea-charm-quark 
distribution is included in our FORTRAN code for WHIT parametrizations. 

Our parametrizations reproduce the distributions within 5% in the region of 0.001 < 
s < 0.8 (corresponding to the region of 4.01GeV2 < < 207GeV2) and 0.001 <x< 0.99 
for the up, down and gluon distributions. Note that the up and down quark distributions 
are written in terms of the valence- and sea- quark distributions as shown in Eq.( p.5|) . In 
the region < s < 0.001 (4 GeV^ < < 4.01 GeV^), they do not agree precisely with 
the distributions mainly because our over-simplified initial distributions do not satisfy the 
perturbative relations among quarks and gluon distributions such as Eq. (p^) . As for the 
charm quark distributions, our convolution integral routine is sufficiently accurate that their 
error is at most 5% reflecting the error in the parametrizations of the gluon distributions. 

2. Parametrizations in the region lOOGeV^ < < 2500GeV2 

In this region, we parametrize the solutions of the massive inhomogeneous AP equa- 
tions obtained in Sec. 0, in which the charm quark mass is retained. According to the 
prescription of Sec. we give the deviations from the QPM predictions for the valence- 
charm-quark and the sea-charm-quark distributions that are defined in Eq. ([4.7|) , in addition 
to the valence-light-quark, the gluon and the sea-light-quark distributions. 

The valence-light-quark, the gluon and the sea-light-quark distributions are parametrized 
in exactly the same manner as in the previous lower case, except that Qq and A in the 
definition of s are now changed to Ql = lOOGeV^ and A = A5 = 302.3 MeV. 

The deviation in the valence-charm-quark distribution is parametrized as 

x6c^{x,s)/a = {Ao + Aix + Agx^ + A3X^)x^{l - x)'^ . (B6) 

Note that, as seen from Eq. ( |4.8a| ), Sc^ is completely universal, i.e. it is common for all six 
WHIT parton distributions. The parametrization of the deviation in the sea-charm-quark 
distribution takes the following form: 
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x6csea{x,s)/a = A x^^'>+^'''\l - xf, (B7) 

which is the same form as the parametrization of the sea-hght-quark distribution. The total 
charm quark distribution is given by Eqs. ( |4.5| ), ( ^.Tp . The QPM part of the sea-charm- 
quark distribution is calculated by using the parametrization of the corresponding gluon 
distribution. 

The relative error of these parametrizations is less than 5% for the up, down and gluon 
distributions in the region of < s < 0.4 (corresponding to the region of lOOGeV^ < < 
3120GeV^), and 0.00125 < a; < 0.99. Note that the up and down quark distributions are 
related to the valence- and the sea- light-quark distributions by Eq. (|2.5| ). For the charm 
quark distributions, the accuracy is also within 5% in the same region but in the different 
X region: 0.00125 < x < 0.99/a, where x=l/a = l/(l + Aml/Q'^) represents the threshold 
for the charm quark. 

3. Prescription in the region of < 4GeV^ 

Finally, we give a prescription for the parton distributions at < Qq = 4 GeV^. This 
is because occasionally one wants to estimate the effects at lower region. We give a very 
crude estimate here that can be used in such cases, rather than setting all the distributions 
to zero or freezing the scale-dependences. 

As is explained in the main text, we give initial parton distributions at Q"^ = Qo = 4GeV^ 
and evolve them to higher by using the inhomogeneous AP equations. Since the AP 
equations cannot generally be solved in the backward direction due to its instability, there 
are no ways to calculate the correct parton distributions below 4GeV^. 

Our prescription to estimate the parton distributions in the region of < 4GeV^ 
is simply to multiply the factor log(Q^/A|)/ log(4GeV^/A4) to the corresponding parton 
distributions at = 4GeV^ i.e. the input parton distributions. This crude prescription 
gives a reasonably good estimate for the valence quark distributions, but the resulting gluon 
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distribution may not be realistic. Therefore, one has to be careful when using this estimate 
at < 4 GeV^. 

In our application, the charm quark pair production process probes the photon structure 
down to ml, which can take the lowest value of (1.3GeV)^ = 1.69 GeV^ for rUc — l.SGeV. 
We have checked that the effect of modifying the prescription at < 4 GeV^ is neghgibly 
small for the total charm quark production cross section. 
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TABLES 



TABLE L The data of F2 adopted in the fit of the valence-quark parameters. 



colUder 



collab. 



(Q2)(GeV2) 



X bins 



Ref. 



PETRA 



PEP 



PLUTO 



TASSO 
JADE 

TPC/27 



4.3 0.03-0.17, 0.17-0.44, 0.44-0.80 flj] 

9.2 0.06-0.23, 0.23-0.54, 0.54-0.90 ^ 

45.0 0.25-0.50, 0.50-0.75, 0.75-0.90 H 

23.0 0.20-0.40, 0.40-0.60, 0.60-0.80, 0.80-0.98 |l| 

24.0 0.10-0.20, 0.20-0.40, 0.40-0.60, 0.60-0.90 H 

100.0 0.10-0.30, 0.30-0.60, 0.60-0.90 |l| 



5.1 0.02-0.20, 0.20-0.36, 0.36-0.74 

20.0 0.196-0.386, 0.386-0.611, 0.611-0.963 



m 

Hi 



TRISTAN AMY 

TOPAZ 



VENUS 



73.0 0.125-0.375, 0.375-0.625, 0.625-0.875 

5.1 0.076-0.20 g 

16.0 0.15-0.33, 0.33-0.78 H 

80.0 0.32-0.59, 0.59-0.98 || 

40.0 0.09-0.27, 0.27-0.45, 0.45-0.63, 0.63-0.81 H 

90.0 0.19-0.37, 0.37-0.55, 0.55-0.73, 0.73-0.91 H 



LEP 



OPAL 



5.9 0.091-0.283, 0.283-0.649 

14.7 0.137-0.324, 0.324-0.522, 0.522-0.836 
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TABLE II. The minimal and the valence-quark parameters as obtained by the best fit. 
Whenever available, we have taken into account correlation in errors. Degree of the freedom of the 



fit is 47-3=44. 



gluon 


best fit 


valence-quark parameti 


ers 




correlations 












C„ 






P(S„,C„) 


0.5 3 


51.6 


0.930(79) 


0.50(17) 


0.24(25) 


-0.52 


-0.75 


0.88 


0.5 6 


53.2 


0.933(78) 


0.51(17) 


0.28(26) 


-0.52 


-0.75 


0.88 


0.5 9 


54.0 


0.938(78) 


0.49(16) 


0.28(25) 


-0.52 


-0.75 


0.88 


0.5 15 


54.3 


0.948(78) 


0.44(16) 


0.26(25) 


-0.52 


-0.75 


0.88 


1.0 3 


54.0 


0.873(76) 


0.77(21) 


0.41(29) 


-0.52 


-0.73 


0.88 


1.0 6 


58.2 


0.882(74) 


0.77(20) 


0.48(30) 


-0.53 


-0.74 


0.89 


1.0 9 


60.2 


0.892(74) 


0.71(19) 


0.47(29) 


-0.53 


-0.74 


0.88 


1.0 15 


60.4 


0.911(75) 


0.61(17) 


0.42(28) 


-0.53 


-0.75 


0.88 


1.5 3 


59.8 


0.821(73) 


1.11(26) 


0.62(34) 


-0.53 


-0.72 


0.89 


1.5 6 


67.8 


0.837(71) 


1.05(24) 


0.70(34) 


-0.54 


-0.73 


0.89 


1.5 9 


70.9 


0.853(71) 


0.95(22) 


0.67(33) 


-0.54 


-0.73 


0.89 


1.5 15 


70.1 


0.879(72) 


0.79(19) 


0.58(31) 


-0.53 


-0.74 


0.89 


TABLE III. 


The values with the standard valence-quark parameters. 




name 


gluon 




name 




gluon 














Ag 


Cg 


X^ 




0.5 


1 


50.7 




1.0 


1 


54.2 




0.5 


2 


51.1 




1.0 


2 


55.0 


WHITl 


0.5 


3 


51.6 


WHIT4 


1.0 


3 


56.0 




0.5 


4 


52.2 




1.0 


4 


57.0 




0.5 


5 


52.7 




1.0 


5 


57.9 




0.5 


6 


53.2 




1.0 


6 


58.6 




0.5 


7 


53.6 




1.0 


7 


59.3 




0.5 


8 


54.0 




1.0 


8 


59.8 


WHIT2 


0.5 


9 


54.3 


WHIT5 


1.0 


9 


60.2 




0.5 


10 


54.5 




1.0 


10 


60.6 




0.5 


11 


54.7 




1.0 


11 


60.9 




0.5 


12 


54.9 




1.0 


12 


61.1 




0.5 


13 


55.1 




1.0 


13 


61.2 




0.5 


14 


55.2 




1.0 


14 


61.3 


WHITS 


0.5 


15 


55.3 


WHIT6 


1.0 


15 


61.4 




0.5 


16 


55.3 




1.0 


16 


61.5 
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TABLE IV. CoefHcients of the parametrization for WHITl parton distribution in the photon. 



Q2 4GeV2 < Q2 < lOOGeV^ 










1 


2 




4 


1v 


Ao 


1.882 


1.213 


0.697 










Ai 





—2.361 


-1.136 










A2 





0.528 


2.406 










B 


0.500 


0.02107 


0.00413 










C 


0.2500 


—0.2376 


0.2018 


—0.0504 





9 


A 


2.000 


-3.28 


2.894 


-1.561 


0.818 




B 





-0.761 


-0.0490 


0.446 







C 


3.000 


1.586 


-0.949 


2.425 












0.461 


0.1041 


—0.01753 


—0.2717 




A'l 





0.00968 


—0.417 


-0.395 


0.843 




B' 


—0.414 


—0.0606 


0.2847 


-0.507 












U.ooo 


1 00c 


U.e5Uy 


U 


Qsea 




0.651 


1.291 


—4.47 


5.14 


— 2.091 




-DO 




0.0901 


—1.356 


1.582 


—0.644 






2.084 


7.74 


—29.70 


38.6 


—17.05 




c 


7.00 


— 16.08 


46.7 


—57.1 


23.86 




Q 






1 f\f\/~1 IT'S ^ ^ OPTtn/"! T 7'2 

lOOCjteV < Q < 2500GeV 












1 


2 




4 


Qv 


Ao 


3.058 


2.474 


1.002 










Ai 


—2.182 


—4.48 


—0.2251 










A2 


1.522 


4.31 


1.314 










B 


0.517 


0.0404 


—0.02100 










C 


0.1655 


—0.02062 


0.0536 








9 


A 


0.784 


—2.238 


16.17 


—62.5 


83.9 




B 


-0.403 


-1.307 


8.78 


-35.8 


53.5 




C 


4.45 


1.027 


44.6 


-160.0 


181.6 






0.3010 


1.275 


— 1.563 


4.10 


-13.37 






—0.1305 


— 1.245 


2.438 


-2.539 


12.73 




B 


—0.489 


0.955 


—4.40 


10.22 


—17.13 




c 


1.331 


—0.2481 


1.950 


—2.072 





Qsea 


A 


0.625 


—0.589 


4.18 


—12.06 


12.57 




Bo 


—0.2492 


—0.411 


0.966 


—2.584 


2.670 




Bi 


2.100 


—5.75 


47.8 


—140.7 


147.6 




C 


4.78 


4.86 


—48.9 


147.7 


—160.2 




Ao 





0.1219 


6.20 


-25.04 


30.98 




Ai 





1.913 


-76.9 


318. 


-392. 




A2 





—7.16 


250.3 


-1062. 


1308. 




A3 





3.19 


—230.1 


1012. 


—1250. 




B 


0.499 


3.47 


—15.26 


19.67 







C 


0.329 


8.24 


-38.0 


46.3 







A 





-0.01815 


0.002043 


-0.00413 







Bo 


-0.3086 


-0.2565 


0.0984 










Bi 


1.376 


-0.463 


1.232 










C 


3.65 


0.729 


-7.57 


7.79 






TABLE V. CoefHcients of the parametrization for WHIT2 parton distribution in the photon. 



Q2 4GeV2 < Q2 < loOGeV^ 







s° 


5! 


s2 


s3 






Ao 


1.882 


1.213 


0.697 










Ai 





-2.361 


-1.136 










A2 





0.528 


2.406 










B 


0.500 


0.02107 


0.00413 










C 


0.2500 


-0.2376 


0.2018 


-0.0504 





9 


A 


5.00 


-14.99 


26.17 


-25.30 


10.12 




B 





-0.937 


0.410 


0.0339 







C 


9.00 


0.709 


3.118 


-0.000582 












0. Kil 


0.1011 


-O.Oirr):-! 


-0.2717 









0.00968 


-0.417 


-0.395 


0.843 




B' 


-0.414 


-0.0606 


0.2847 


-0.507 







C 


1.244 


0.588 


-1.228 


0.809 





Qsea 


A 


1.237 


3.39 


-10.75 


12.46 


-5.58 




Bo 


—0.0727 


0.1748 


— 1.392 


1.711 


—0.796 




D 


A on 


17.87 


—58.1 


81.9 


At A 

— 41.4 






±4.o4 


— 44. y 


iiy. / 


1 K 

— locS.o 


/ 0.0 










lOOGeV^ < < 2500GeV2 










s° 




s2 


s3 




Qv 


Ao 


3.058 


2.474 


1.002 










Ai 


-2.182 


-4.48 


-0.2259 










A2 


1.522 


4.30 


1.315 










B 


0.517 


0.0403 


-0.02098 










C 


0.1655 


-0.02063 


0.0537 








9 


A 


1.095 


-2.388 


9.19 


-30.32 


34.8 




B 


-0.441 


-0.907 


4.68 


-18.66 


27.17 




C 


10.99 


4.71 


28.01 


-127.9 


164.0 






0.3010 


1.275 


-1.563 


4.10 


-13.37 




^'1 


-0.1305 


-1.245 


2.438 


-2.539 


12.73 




B' 


-0.489 


0.955 


-4.40 


10.22 


-17.13 




c 


1.331 


-0.2481 


1.950 


-2.072 





Qsea 


A 


1.287 


-2.069 


11.57 


-35.7 


37.4 




Bo 


-0.2340 


-0.443 


1.235 


-3.72 


3.84 




Bi 


6.46 


-10.48 


89.8 


-284.7 


299.8 




(• 


'■,.?>:> 


10.11 


-l.'i'l.V 


127. 


- LIV. 


SCv 


Ao 





0.1219 


6.20 


-25.04 


30.98 




Ai 





1.913 


-76.9 


318. 


-392. 




A2 





-7.16 


250.3 


-1062. 


1308. 




A3 





3.19 


-230.1 


1012. 


-1250. 




B 
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TABLE VI. CoefHcients of the parametrization for WHITS parton distribution in the photon. 
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TABLE VII. Coefficients of the parametrization for WHIT4 parton distribution in the photon. 
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TABLE VIII. Coefficients of the parametrization for WHIT5 parton distribution in 


the photon. 
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TABLE IX. Coefficients of the parametrization for WHIT6 parton distribution in 


the photon. 


V 


















i> 




«2 

o 


o 


s 


Qv 


^0 




9 nnn 

Z.UUU 


n 71 Q 
U. ( 10 


n 
U 









U.UDzo 


7 ni 


n 1 9c^i 
U. IzOl 


n 
U 







A^ 




— U.400 


in AQ 

lU.4cS 


9n 











n n97Qfi 
— u.uzr / yo 


— U.UUOUu 


n 
u 







n 




— ±.zoo 


1 .y^:! 


— u.yyo 





9 


A 

yi 


io.UU 


— Di.U 


1 97 fi 


1 Q 

— loy.y 


59.9 




D 

ID 


U 


1 1 HQ 


n Q/1 K 
U.o40 


— U.oOl 











u.ioyo 


4. lo 


n 1 7Kc; 
— U.l /DO 









() 


0.()02 


()..)()!) 


— 2.0;> i 


l.o!)2 




A' 


n 
U 


n nQ99 
— u.uyzz 


1 SQQ 

— i.csyy 


^.lo 


—2.494 




B' 




U.O 1 u 


— 1 71 Q 
— 1.1 ly 


1 116 







O' 




— U.OO 1 


u.oou 


U. ( oo 


— U. lOZ 


Qsea 


A 
J\ 


o. io 


o.oy 


99 C7 


1 Ofi 

lo.yD 


—5.14 




Rex 


—0.1003 


0.1603 


— 1.037 


0.944 


—0.2915 




Bi 


5.69 


18.67 


-46.7 


50.5 


— lO.oO 




C 


21.49 


-56.5 


129.3 


-145.9 


57.5 










innnpv^ < < 9fsnnopv2 

XUUVjrt; V \ li^ \ ZUUUVjrt; V 










„0 
o 


«ji 


^2 
o 


b 


s 


Qv 


A^ 


41:. z ( 


o.uyo 


1 K91 


n 
U 







A , 


A '7 A 
— 41:. / 1 


on 


9 /tQQ 


U 







J±2 


z.oo < 


Afi 


/I 1 n 

4. iU 


U 









U.D / 


n O'iQzl 
— U.Uoy4 


u.ui 1 oo 


n 
U 









n 1 709. 


— u-uZr-yryo 


U. 14:01 


n 

u 





9 


A 
>1 


Z.OiO 




U.OoO 


o.o4 


—9.92 




D 


— U.4 ( y 


n fin7 
— U.uU / 


1.400 


nQ 


9.33 






i ( .uo 


41:. yo 


OA Q7 

z4.y ( 


1 C^Q 9 


295.4 




4' 
^0 




1 1 7Q 

i.i 1 y 


1 Q1 

— i.yio 


7 9^; 
1 .zo 


—18.39 




^1 


— U.lDUU 


1 11/1 

1. 1141: 


9 QQQ 


— U.DO 


19.23 








U. 1 oo 


— O.oU 


1 n 7'\ 

lU. 1 o 


—19.93 






1 9Q7 


n 1 (^ftQ 
— u. iDoy 


1 .yuD 


9 01^7 
— Z.UO ( 





Qsea 




o.o4 


c; fil 
— U.Dl 


OU.U 


oon 7 


302.8 








— u.4uy 


9 Ofi'i 


1 n c;n 
— lU.OU 


14.87 






7Q 


— o.oo 


1 fiA n 


71 9 

— / iz. 


973. 




y 


J. i{) 




— Z i o . I 


1 1 7", 


— 1.)!)2. 


OCv 


^0 


n 
U 


n 1 91 Q 
u.iziy 


A 9n 


9c; nA 


30.98 






n 
U 


1 01 Q 

i.yio 


7fi O 

— ( D.y 


Q1 Q 

olO. 


—392. 




A^ 


n 
U 


7 1 

— / .ID 


9c;n "i 


1 nfi9 
— lUDz. 


1308. 




A ^ 


n 
U 


Q 1 O 

o. iy 


9Qn 1 


1 ni 9 

lUlz. 


—1250. 




D 


n zlQQ 


'i A7 
0.4 ( 


1 Pi 9(^ 


1 Q fi7 

ly.o ( 









u.ozy 


8 24 


QQ n 


Afi 


n 
U 


ocsea 


/I 


U 


n nAc\n 

— u.U4yy 


n 1 nofi 


n n7Q7 







Bo 


-o.;^tii 


-().r,70 


0.22o7 












7.68 


-8.83 


3.88 










r 


2.r,is 


O.GOl 


-8.70 


10. Go 


n 



FIGURES 

FIG. 1. The data and the theoretical predictions for the photon structure function F2 /a. The 
vertical axes show F2 /a and the horizontal axes represent the scaling variable x. The data points 



with a cross mark are not used for the fit, according to the selection criterion Eq. (3J). The OPAL 
data points are obtained by removing the direct charm quark contribution. We hence drop from 
the theoretical curves the 7*7 cc contributions. The remaining structure function data contain 
all hadronic final states. (a) The best fits with WHIT1~WHIT3 gluon distributions, (b) The 
best fits with WHIT4-WHIT6 gluon distributions. 

FIG. 2. The deviation of each data point from the best fit, -^2^(2;, Q^)fit — -^2*^(2;, Q^)data divided 
by the error o"(F2^(x, (5^)data)- A data point with smaller x value is placed lower in each data set 
with a common (Q^). The data points with a simple cross mark are removed from the fit by the 
criterion Eq. (3.1). 



FIG. 3. The Cg dependence of the total for fixed Ag^s when the standard valence-quark 
distributions are taken for Ag = 0.5 (solid line) and 1.0 (dashed line). The large square, diamond 
and cross marks are the minimal values as obtained by tuning the valence-quark distributions. 

FIG. 4. Gluon distributions at = 4^ 20 and 100 GeV^. The top 3 figures are for WHITl 
(solid), WHIT2 (dotted) and WHIT3 (dashed), the middle 3 figures are WHIT4 (sohd), WHITS 
(dotted) and WHIT6 (dashed), and the bottom 3 figures are GRV [|] (solid), DG @ (dotted) and 
LACl [|l^ (dashed) for comparison. 

FIG. 5. Predictions for the charm quark distribution in the photon, calculated by the QPM 
(solid lines) and by the massive inhomogeneous AP equations (dash-dotted lines). The va- 
lence-up-quark distributions of WHIT1-WHIT3 (dashed lines) and WHIT4-WHIT6 (dotted hues) 
are also shown for comparison. 



FIG. 6. The differential cross section of the process e^e" e'^e'cc at TRISTAN en- 
ergy {y/s = 58GeV), evaluated exactly and via EPA: (a) near the charm quark pair threshold 
= rricc = 4GeV) and; (b) at far above the threshold = rricc = 15GeV). Vertical bars 
indicate errors of numerical integration of the exact matrix elements. We set mc=1.5 GeV and 
a=l/137. 

FIG. 7. The invariant mass distribution of charm quark pair in the process e"'"e~ e^e~cc 
evaluated exactly and via EPA at -s/s = 58GeV. We set rric = 1.5 GeV, and a = 1/137 

FIG. 8. The leading order prediction for the total cross section of the inclusive process 
e+e" e^e~ccX. The contributions from the resolved photon processes depend on the par- 
ton distributions WHITl to 6. The curves are obtained by setting rric = 1.5GeV, a = 1/137 and 
A4 = 0.4 GeV for the strong coupling, and by requiring the invariant mass W of the hadron system 
to satisfy W>2mD = 3.74 GeV. The vertical bars attached to the WHITl, WHIT4 and WHIT6 
predictions indicate the dependence of the cross sections on the charm mass rric between 1.3 GeV 
and 1.7 GeV. 

FIG. 9. The leading order predictions for the total cross sections of the inclusive process 
e^e" — ^ e+e^ccX, where contributions of the direct and resolved photon processes are shown 
separately. The parton distributions of WHITl (a) and WHIT6 (b) are used to calculate the 
resolved photon contributions. We set rric = l.SGeV, a = 1/137 and A4 = 0.4 GeV. The vertical 
bars indicate the dependence of the cross sections on the charm mass rric between 1.3 GeV and 1.7 
GeV. 
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